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Abstract 

We study the low-bias stabilization of active magnetic bearings (AMBs) subject to voltage saturation based on a 
recently proposed model for the AMB switching mode of operation. Using a forwarding-like approach, we construct 
a stabilizing controller of arbitrarily small amplitude and a control-Lyapunov function for the AMB dynamics. We 
illustrate our construction using a numerical example. 

1 Introduction 

Active magnetic bearings (AMBs) are being employed in a variety of rotating machines (e.g., compressors, milling 
spindles, and flywheels) in place of conventional mechanical bearings. Due to the non-contact nature of the magnetic 
bearings and rotor, AMBs have the unique ability to suspend loads with no friction, operate rotors at higher speeds, and 
operate in applications where the use of lubricants is prohibited. Since AMBs can be actively controlled, they offer 
other potential advantages over mechanical bearings, viz., the elimination of vibration through active damping, the 
adjustment of the stiffness of the suspended load, and the automatic balancing of rotors. AMBs are normally operated 
with all electromagnets active at all times. An alternative mode of operation is to activate only one electromagnet 
along each direction at any given time. Although posing more difficulty to the control design, this switching operation 
prevents opposing electromagnets from producing counteracting forces, and thus helps reduce power consumption. 

Typically, an AMB is operated by introducing a sufficiently high, fixed magnetic flux in each electromagnet, which 
is referred to as the bias flux. The bias value is normally set to a fraction of the saturation flux of the electromagnet. 
This procedure facilitates the design of the AMB 'control' flux, which is superimposed on the bias flux. Specifically, 
this conservative practice allows the system to be modeled by a controllable linear system, thus, enabling the use of 
linear control design techniques; see for example G][7l[9). Although the bias flux facilitates the control synthesis, it 
increases electric power losses in the AMB system, causing rotor heating and affecting the machine efficiency. While 
lowering or eliminating the bias flux is desirable in order to minimize power losses, it enhances the AMB system 
nonlinearities and may lead to a control saturation 1 or singularity. Due to these conflicting objectives, the design of 
AMB controllers with reduced power loss is a challenging problem. 

In this paper, we consider the problem of low-bias 2 control of AMBs operating in the switching mode with the 
constraint that the input voltages are amplitude limited. This problem was previously addressed in 1151 1161 using 
the nested saturation design method of Teel 1131 II 41 . In addition, an optimal solution was sought to the power-loss 
minimization problem under voltage saturation in |2|. (A comprehensive literature review of low-, asymptotic-zero- 
, and zero-bias AMB controllers without voltage saturation can be found in I10II17I .) Here, we pursue a different 
approach to stabilizing the AMB model of 1 16 1. Namely, we use a forwarding-like method 1 8 1 to design a control law 
of arbitrarily small amplitude that renders the AMB system globally asymptotically stable (GAS) to the origin. The 
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main differences between our approach and 1 16 1 are: (i) we allow the use of different saturation functions (both 'hard' 
and 'soft'), thus providing more flexibility in the control implementation, and (ii) we provide an explicit construction 
for a 'global' control Lyapunov function (CLF) for the system. This contrasts with the result of 1161 . which relies on 
the standard 'hard' saturation and whose Lyapunov-like function is only an 'asymptotic' CLF for the system, i.e., a 
CLF only when the system operates in a certain region of the state space. 

The rest of this paper is organized as follows. In Section [2] we introduce the AMB model and motivate our 
stabilization problem. In Section [3] we prove a lemma that constructs a stabilizing controller of arbitrarily small 
amplitude and a CLF for a three-dimensional chain of integrators. In Section|4] we use our lemma to construct a CLF 
and corresponding stabilizing feedback for the AMB model. We provide a numerical example in Section |5] and we 
close in Section|5]with a summary of our work. 



2 AMB Model and Problem Statement 

The original nonlinear electromechanical model of the one degree-of-freedom AMB system shown in Figure 1 below 
can be subdivided into the mechanical subsystem dynamics, the magnetic force equation, and the electrical subsystem 
dynamics. The mechanical subsystem is governed by 



my 



where m is the rotor mass, j/gR represents the position of the rotor center, 6 R is the magnetic flux in the ith 
electromagnet, F{ el denotes the force produced by the ith electromagnet, given by [18) 

Fi = ^ 1 . \ 1 = 1,2, (2) 

Ho is the permeability of air, and A is the cross-sectional area of the electromagnet. The electrical subsystem is 
governed by the equations 1 1 8 1 

NQi + RiU =ih, t = 1,2, (3) 

where N denotes the number of coil turns in the electromagnet, R4 is the resistance of the ith electromagnet coil, 
Vi G K is the input control voltage of the ith electromagnet, Ii e R is the current in the ith electromagnet which is 
related to the flux according to [ 18 1 



2 (flb + (-!)*» )*i 
h = — ; , i = 1,2, (4) 

and go is the nominal air gap. 

In this paper, we use the new form proposed in [ 16 17 1 for the AMB dynamics in the switching mode of operation. 
We only outline the model derivation here, and refer the reader to 1171 for the details and justification. Consider that 
the flux is given by 

$i = $ + &, i = l,2 (5) 



where $0 > is the constant bias flux and <pi is the control flux. Let the generalized control flux be defined as 



and consider the voltage switching strategy 



(6) 



(7) 



Vi = v, v 2 = if 4> > 
Vi =0, V2 = —v if <f> < 

where v is the generalized control voltage. Based on the AMB model ([O-Q has the equivalent form 

y=^A^ + ^\), (8) 
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where $o = + min {</>i (0), 02(0)} and the coil resistance was neglected for simplicity. 

Now, consider that the input voltages to the original AMB model are amplitude limited, i.e., < v m ax, i = 1,2 
where v max is the known limit. Defining the states X\ = y, x 2 = y, and x 3 = <fi and the change of input v = Nu, the 
above model becomes 

±1 = x 2 , x 2 = f3 x 3 + fcxa \x 3 \ , x 3 = u, (9) 

where (3 — 2<&o/(m/i ^4) and j3\ — l/(m/i A), with the input constraint \u\ < v max /N. 

Our goal in this work is to design a feedback control law u(x), satisfying the above saturation constraint, such 
that x — is GAS where x = (xi, x 2 , x 3 ) T . The GAS property means there is a continuously differentiable (C 1 ) 
function V : K 3 — > [0, 00) that is radially unbounded and zero only at the origin, and for which the derivative V along 
all trajectories of the system in closed-loop with the controller u(x) is negative definite. A function V satisfying these 
requirements for some feedback u(x) is called a CLF for l|9} [5 1. See 1 17 1 for a discussion on the relation between the 
global asymptotic stability of x and the stability of the original AMB states. 



3 Preliminary Result 

In the following lemma, we use a forwarding-like approach based on 1 8 1 to construct a CLF and corresponding stabi- 
lizing state feedback of arbitrarily small amplitude for a three-dimensional chain of integrators. We let a denote the 
standard saturation projecting K onto [—1, +1], i.e., u(s) = s when \s\ < 1 and a(s) = sign(s) otherwise. However, 
see Remark 1 below for a version of our construction where the hard saturation a is replaced by a smoother function. 

Lemma 1 Consider the chain of integrators 

6=6, 6=6, 6 = /i, (10) 

where e lis the control input. Let C\ , c 2 > be given and 

S < min<j -,ci }> . (11) 
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Define the change of variables 



and the functions 



zi=S (ci6 + 6) H — 6, z 2 = ci£ 2 + 6, z 3 = 6 (12) 

Cl 



U(z 2 ,z 3 ) = \zl+ [ a(s)ds (13) 
z Jo 

g(zi,z 2 ,Z3) = — cr{zi) +4(1 + U(z 2 ,z 3 )) 

Cl 

x (z 3 + a(z 2 )) - - (z 2 + z 3 ) . (14) 



Then the function 



^(6,6,6) = W(z 2 ,z 3 ) + 2U 2 (z 2 ,z 3 ) 



s 



-■^z 2 z 3 + J a(s)ds (15) 

is a CLF for dlOl since its derivative along the trajectories of system dlOl l in closed loop with the control 

^(6,6,6) = -c\a{z 2 ) - c 2 a(g{zi,z 2 ,z 3 )) (16) 
is negative definite. In particular, the bounded control (TTBT l renders ( fTol GAS to the origin. 
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Remark 2 The saturation function a from Lemma[2is Lipschitz continuous but not continuously differentiable. We 
chose this hard saturation for simplicity of analysis. However, one can prove a similar result with smoother saturation 
functions. For example, if one instead chooses the C 1 function 



a(s) 



s if \s\ < 1 

sign( S ) [1 + i tan" 1 (vr (|s| - 1))] if |s| > 1, 



then the only change in the statement of Lemma^is to the right-hand side of the inequality il It . The proof for this 
saturation is similar to the proof we give below. 

Remark 3 Condition il y in Lemma^allows cases where c\ is arbitrarily large and 5 is as small as desired, as well 
as cases where max{ci, c 2 , 5} is arbitrarily small, including the case where c\ = c 2 = S G (0, 4/3]. 

Proof. We use d!2t and /i = — c\a{z 2 ) + a, where a is the new control input, to transform i ll Ob into 

zi = 5 (z 2 - <j(z 2 )) + -^a 

zi = c\z 3 - aa(z 2 ) + a (17) 
z 3 = — Cia(z2) + a. 

We first perform a Lyapunov-type analysis on the (z 2 , ^-subsystem using the function 

V 1 (z 2 ,z 3 ) = 4U(z 2 ,z 3 ) + 2U 2 (z 2 , z 3 ) - -z 2 z 3l (18) 

where U(z 2 , z 3 ) was defined in (1131 . By separately considering the cases where \z 2 \ > 1 and \z 2 \ < 1, one can easily 
verify that Jl 81 is positive-definite and radially unbounded for S < 2. In fact, if \z 2 \ < 1, then 

Vi(z 2 ,2 3 ) > 4U(z 2 ,z 3 )-\z 2 z 3 \ > 2z 2 + 2z 2 - \z 2 z 3 \ 

> \{4 + 4)- d9) 

On the other hand, if \z 2 \ > 1, then J* 2 a(s)ds = 1/2 + (\z 2 \ - 1) > |z 2 | /2, so 

^1(22,^3) > 4 (^i + I l^l) +l{^4 + \\z2 

- \Z2Z 3 \ 

> 2z 2 + 2\z 2 \ + -z 2 -\z 2 z 3 \>-z 2 

+2\z 2 \. (20) 

The time derivative of (1181 along the trajectories of the (z 2 , ^-subsystem is given by 

Vx = 4(l + U(z 2 ,z 3 ))(- Cl cr 2 (z 2 ) + (z 3 + a(z 2 ))a) 
-- (ciz 3 - cid{z 2 ) + a) z 3 

-|^2 {-cia{zi) + a) 

2, \ Sci 2 



-4 Cl (l + U(z 2 ,z 3 ))a'(z 2 )-^-z£ 

da , Set 
+—a(z 2 )z 3 + —z 2 a(z 2 ) 

+ [4(l + U(z 2 ,z 3 ))(z 3 + cr(z 2 )) 



-—z 3 z 2 

2 2 



(21) 
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Since \<j{z 2 )\ < \z 2 \ and \z 2 a(z 2 )\ = z 2 a(z 2 ) give 



\ z l + \^ 2 )z 3 < ^CT 2 (z 2 ) < -Z 2 a{z 2 ), 



(22) 



d2Ti becomes 



V, < 



-4 Cl (l + U(z 2 ,z 3 ))a 2 (z 2 )-^-z 2 
+ ^5c 1 z 2 a(z 2 ) 

+ [4(l + U{z 2 ,z 3 )){z 3 + a(z 2 )) 

s s ■ 

— z 3 z 2 

2 J 2 

We now distinguish between two cases: 

1. |z 2 < 1. Then a(z 2 ) = z 2 , z 2 a(z 2 ) = a 2 (z 2 ), and since S < 4/3, we get 



(23) 



Vi < -3 Cl {l + U(z 2 ,z 3 ))a 2 (z 2 ) 



Sci 



+ [A(l + U(z 2 ,z 3 ))(z 3 + a{z 2 )) 



-2 Z3 ~2 Z2 



(24) 



2. \z 2 \ > 1. Then ct(z 2 ) = sign(z 2 ), z 2 cr(z 2 ) = M, and 



(l + [/(z 2 ,z 3 ))a 2 (z 2 ) > 1 + f* 2 a(s)ds 



> i + sds + £ 



l + / sisn(z2) S d S + sign(z 2 )/ f 



z 2 

sign(z 2 ) 



<j(s)ds 



sign(z 2 ) 



fi.S 



(25) 



^2 



Since (5 < 4/3, applying d25t to ( I23l l again yields &241 . Thus ( I24t holds for all 



Z 2 . 



Finally, we perform a Lyapunov-type analysis on the whole system dl7l > using the positive definite radially un- 
bounded function V from dl5l >. Since V{z\, z 2 , z 3 ) — V\(z 2 , z 3 ) + J* 1 cr(s)ds, where V\(z 2 , z 3 ) was defined in 1.181 . 
it follows from 1241 that the time derivative of dl5> along the system trajectories satisfies 



V < -3 Cl (l + U(z 2 ,z 3 ))a 2 (z 2 )-^-zl 

+S\z 2 - <j{z 2 )\ 
5 



a(zi) + 4 (1 + U(z 2 ,z 3 )) {z 3 + a{z 2 )) 



ci 



- — z 3 z 2 

2 2 



(26) 



We again distinguish between two cases: 

1. \z 2 \ < 1. Then z 2 — <r(z 2 ) = 0, so (I26> gives 



V < -2 Ci (1 + U(z 2 ,z 3 ))<j 2 (z 2 )--±z 2 



-a( Zl )+A(l + U(z2,zs)) 



ci 



x (z 3 + ct(z 2 )) - ^z 3 - ^z 2 



(27) 
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2. 1 2a | > 1. Then i25i holds and since \z 2 — o"(^2)| < l 2 2| and 5 < c\, we again get ( 127k 
Thus, ( I27> holds for all z 2 . Using a = —c 2 a(g(zi, z 2) 23)) and dl4> . |27i becomes 

\> < -2c 1 (l + C/(z 2 ,z 3 ))ff 2 (^)-^ 1 ^ 

-C2(r(g(zi, z 2 , z 3 ))g(zi, z 2 , z 3 ) 

< -2 Cl a 2 (z 2 )-^-z 2 

-c 2 a 2 (g(z 1 ,z 2 ,z 3 )) (28) 

by separately considering the cases \g\ > 1 and \g\ < 1, i.e., V is negative definite. One can easily verify that the 
origin is a unique equilibrium point of dl7> in closed loop with a = —c 2 a(g(zi, z 2 , z 3 )). Thus, (z\, z 2l z 3 ) = is 
GAS 1 3 1 . It follows from fl3 that (£1 , £2, £3) = is GAS, and that O is a CLF for 



4 Main Result 

We are now ready to state our main result on the global asymptotic stabilization of (|9j with control saturation. 
Theorem 4 The control law 

»(*) = fl'gf 3 ? . (29) 

where /i(-) was defined in (fT6l and 

-^3 =A)«3+A^3|«3|, (30) 

in closed loop with (|9} ensures a; = is GAS, and has arbitrarily small amplitude. Moreover, with the choice i!5i . the 
function x 1— > V(xi, x 2 , X3) is a CLF for (|9}. 

Proof. We use the change of variable ( I30l > and the change of input ( I29l > to transform (|9jl into the chain of integrators 

£1 = 2:2, i 2 =^3, X s =fi. (31) 

After setting /j, to d!6i with ^1 = x\, ^2 = ^2, and £ 3 = X3, we can invoke Lemma^o conclude (x\, x 2 , X3) = 
is GAS for ( I3U . It then follows from J30i that x = is GAS for the system l|9} in closed loop with the feedback 
( I29M30> . The fact that x 1— » V{x\ , x 2 , X3) is a CLF for the AMB dynamics follows from the proof of Lemma^ 



u{x{m < mm < Eipi w >o. (32) 



Remark 5 Notice that the control law j29t has arbitrarily small amplitude since along the closed-loop trajectories 

0o 



Thus, one can choose c\ , c 2 such that 

c\ + c 2 < — — — (33) 

to ensure the AMB voltages satisfy the voltage constraint \vi(t)\ < v max , i = 1, 2 Vt > 0. From (1331 . one can also see 
that the control gains can be adjusted to accommodate sufficiently small bias levels due to the direct dependency of (3q 
on the bias flux $0 (i.e., a low-bias control). However, one cannot set $0 = since the control may have a singularity. 

Remark 6 It is worth comparing our construction from Theorem|4]with the AMB controller from Section VI of 1 16 1, 
where the stabilizing controller took the form of nested hard saturations. In 1 16 1, the proof that the control stabilizes 
the AMB system is based on the method of 1 13 1, which uses the hardness of the saturation o in an essential way. 
The argument shows that the nested hard saturations cause the AMB dynamics l|9) to assume a simplified form after 
sufficiently large time. Then one shows that this new system of 'asymptotic' equations is GAS using a Lyapunov-like 
analysis. Unfortunately, the Lyapunov function used in the analysis of 1 16 1 is not a CLF for the AMB dynamics (|9jl. 
On the other hand, our construction is not restricted to the standard hard saturation as explained in Remark 1, and thus 
can lead to smoother stabilizing controllers. Also, the Lyapunov function from Theorem|4]is negative definite along 
the closed-loop AMB trajectories with the feedback d29t for all time, and thus is a CLF for the AMB dynamics l|9}. 
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Remark 7 Knowing a global CLF V(x) and corresponding globally stabilizing feedback u(x) for the AMB dynamics, 
we can construct a feedback K{x) that renders the system input-to-state stable (ISS) to exogenous disturbances II II . 
Choosing V(x) and u(x) from our theorem, and letting x = F(x, u) denote the open-loop AMB dynamics, it follows 
from [9, Theorem 1] that the system x = F(x, K(x)+d) is ISS relative to relative to any measurable locally essentially 
bounded (e.g., locally bounded piecewise continuous), exogenous disturbance d(t) under the combined feedback 
K(x) = u(x) — dV(x)/dxz. This is because the term dV(x)/dxs in the feedback is the Lie derivative of V in 
the direction of the constant vector field x i — > (0, 0, 1) T that multiplies the input in the AMB dynamics. Note that 
in this case the feedback K{x) is not necessarily saturated. This is because if K{x) were any saturated feedback 
rendering x = F(x, d) ISS, and if K were a positive global bound on K(x), then the third component of the dynamic 
±3 = K(x) + d would have unbounded trajectories when the disturbance is for example d(t) = 2K , thus violating 
the ISS requirement. Hence, there is no globally saturated feedback that renders the AMB dynamics ISS to exogenous 
disturbances. Moreover, the CLF allows us to construct a feedback K(x) rendering the noisy system x = F(x, K(x + 
e) + d) ISS with respect to the disturbance d for small measurement errors e(t) in the controller 1121 . In the context 
of AMBs, d could represent rotor unbalance or vibrations transmitted through the bearing foundation while e could 
represent inaccuracies in sensing displacement, velocity, and/or flux. 

5 Numerical Example 

The proposed control law from Theorem|4]above was simulated on the same high-fidelity model of the one degree-of- 
freedom magnetic levitation system 1 4 6 1 used in [ 1 6 1 . The model includes flux leakage, magnetic material saturation, 
coil resistance, voltage saturation, and flexible modes, and thus is a fairly accurate representation of the actual system 
test rig |6|. We refer the reader to |6| for a complete description of the system model. 

As in [16 1, the system was simulated with $o = 10 /iWb, w m ax = 10 V, and initial conditions set to y(0) = 0.15 
mrad, y(0) = mrad/s, $i(0) = 20 fiWb, and $2(0) = 60 /iWb. The above flux values are for the total flux in the 
magnetic levitation system, i.e., including gap flux and flux leakage. The control parameters ci, C2, 6 were selected 
according to (1331 and (II II to yield the fastest convergence possible for the state y. This resulted in the following 
parameter values: c\ = 5465, C2 = 5300, and 5 = 4/3. The simulation results for the system states y, y, and <£> g i 
(i = 1, 2), and for the control voltages v.i (i — 1,2) are shown in Figure 2 below, where $ 9 , denotes the gap flux. 
When compared with the simulation of the controller from Section VI of [16 1 (see Figure 8 of 1 16 1), the our control 
yields similar closed-loop performance although it can be tuned to produce a faster response as seen from Figure 2. 

6 Conclusion 

We constructed a stabilizing controller and control-Lyapunov function for low-bias active magnetic bearings (AMBs) 
with voltage saturation. The AMB dynamics used in the construction is borrowed from a recently-developed model for 
AMBs operating with only one electromagnet active at any given time. Our construction is based on a new forwarding- 
like approach to stabilizing a three-dimensional chain of integrators with input saturation. The proposed AMB control 
law has arbitrarily small amplitude and the flexibility of being implemented with different saturation-like functions. 
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Figure 1 : Schematic of the AMB system. 
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Figure 2: State trajectories and control voltages. 
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